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Abstract 

^^ ' I discuss some compelling suggestions about particles which could be the dark matter in the 

Q\ ■ universe, with special attention to experimental searches for them. 

On 

This is the write-up of a talk given at a conference where I was probably the only particle physicist. My 
task was to give an overview of one area of cosmology which has a symbiotic relation with particle physics 
\& ' — viz., the candidates for dark matter. 

I could proceed in two ways. One is to give very brief outlines of all different suggestions (or all I know 
about) for dark matter candidate particles. The second is to focus on a few compelling ones, and give some 
details of them. After some thought, I chose the second path. I will divide this talk into three sections, by 
three different kinds of particles which are potential candidates for dark matter. 

^ ! 1 Light neutrinos 

/^v \ Neutrinos interact very feebly, only through weak interactions. Despite this fact, they were in thermal 
equilibrium in the very early universe because the density of particles was much higher at that time, and so 
the mean free path (or the average reaction time) was small. At that stage, their number density differed 
from that of the photons only because one is a fermion and the other is a boson, and the relation was 

5— i ' 3 

to! n v = -n~ l . (1.1) 

<3 ■ 

This age of thermal equilibrium ended once the density of the particles became small enough in an expanding 
universe, and it happened when the temperature was about 1 MeV. After that, the neutrinos followed only the 
overall expansion. The photon number density also fell by the cube of the scale factor, so the ratio remained 
the same. The only exception to this statement occurred when the temperature dropped somewhat below 
the electron mass, and e + e~ annihilations produced new photons. This event increased the photon number 
density by a factor of 11/4, so that in the present universe 

Wo^4 1lWo" 110Cm "- (L2) 

If the neutrinos are massive, the energy density is obtained by multiplying this by the mass. The total 
energy density of the universe is parametrized by 10 4 h 2 tt eV/cm , where h is the Hubble parameter in units 
of 100 kms -1 Mpc~ . Thus, the fraction of energy density contributed by light neutrinos would be 

= m v n v = / m v \ 1 

V lO^fteV/cm 3 V92evJ h?Sl ' ( ' j 

Atmospheric neutrino data indicate neutrino masses of order 10" 1 eV. Solar neutrino data indicate even 
smaller masses. Although these do not apply for all neutrino species, it surely is suggestive of the fact that 
F„<1. 
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"Invited talk at the conference "Cosmology: Theory confronts observations", held at IIT Kharagpur, January 1999. 
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Figure 1: Neutrino mass and degeneracy allowed by cosmological energy density bounds. The curves are 
contours for h 2 VtF v = |, \ and ^ respectively. From Ref. ||]. 



Before discarding light neutrinos as dark matter candidates on this ground, I would like to emphasize that 
there is a very strong assumption implicit in this entire argument, viz. that the neutrinos have no chemical 
potential so that their number is the same as that of the antineutrinos. If chemical potential is nonzero, the 
number density of neutrinos depend on it as well as the temperature, and one cannot even start with the 



simple relation of Eq. (1.1). 

If the neutrinos have chemical potential, they can contribute a larger fraction to the energy density of 
the universe. This issue has been analyzed recently Q. We show the results in Fig. |j. The horizontal axis 
corresponds to neutrino mass. The vertical axis is n v — rip in units of the CMBR n 7 in the present universe. 
The lines in this plot correspond to the values h^VLFp = \, \ and \ respectively, starting from the inner 
line. 



2 Axions 

So far we talked about light neutrinos. They are known to exist. The other candidates to be discussed have 
not yet been detected in any experiment. Thus, I need some motivation for talking about them. 

For the case at hand, the motivation comes from strong interactions. These interactions are believed 
to be described by the gauge theory called QCD. In a gauge theory, one decides on some symmetry of the 
Lagrangian and writes down all terms that are cosistent with this symmetry. For strong interactions, this 
symmetry was found in the 1970s. If one cares to do only perturbative calculations with the Lagrangian 
obtained from this symmetry, there is no problem. If, however, one includes non-perturbative effects, one 
encounters one term which violates parity and time-reversal symmetries: 



-QCD 



-pert 



32tt 



r^u^^/j 



(2.1) 



where G^ is the generalization of the electromagnetic field-strength tensor for the relevant symmetry. Here 
g is the gauge coupling constant, and 8 is an arbitrary parameter. 

The term involving this parameter is the generalization of a term E ■ B, which violates parity and 
time-reversal symmetries. Among other things, this term will give rise to an electric dipole moment for 
the neutron. The experimental bounds are very strong on this, which restricts this parameter severely: 
6 < 1(T 10 . 



Very small parameters always present certain conceptual difficulties in quantum field theories. If they 
are not zero, one has to justify them by some symmetry. Peccei and Quinn discovered such a symmetry by 
introducing a Higgs field a(x) which has an interaction 

(j iiv(j jiv , (2-2) 



32tt 2 /p Q 

where /pq has the dimensions of mass. The field, of course, has other interactions, including a potential of 
its own, which is minimum at a — OfpQ- In general, now, we can write a(x) = 0fpQ + a(x). This field a{x) 
vanishes at the minimum and is qualified to be a quantum field. It is called the axion field. 



Notice how it has solved our original problem. If we add the terms shown in Eqs. (2A) and (^2) and use 
the definition of the quantum field, we see that the parameter 9 cancels out. Said another way, 9 is forced 
to equal to zero by this procedure. We have an interaction of the axion obtained by replacing a by a in Eq. 
( |2.2| ), but that does not violate parity and time reversal provided the axion is intrinsically negative under 
these symmetries. 

We now examine how the axion can be important as the dark matter. The equation of motion of the 
axion field, neglecting its interactions, is 

a + 3Ha + m 2 a a = Q, (2.3) 

assuming the field to be spatially homogeneous. This is similar to the Klein-Gordon equation, but there are 
two differences. One of them is the second term on the left side, which shows the effects of the expansion 
of the universe. The other is hidden in the notation, and is the fact that the 'mass' m a appearing in the 
equation is not a constant. It arises due to co-operative effects, and depends on the temperature. 

At early times, when the temperature was very high, m a was zero. This is because all non-perturbative 
effects involve a factor exp(— 1/g 2 ), which is negligible since at high energies g is very small for QCD. In 
this stage, the solution of Eq. (|2.3| ) is given by a — constant. As the universe cools down, g increases, and 
becomes large enough to make the non-perturbative effects important when the temperature drops down 
sufficiently, say below a value Xqcd- To find the solution after this time, let us first forget about H and also 
the time-dependence of m a . In that case, the solution would be 

a = Acosm a t , (2.4) 

where A is a constant. If H <C m a and also the time variation of m a is small, we can still try a solution like 
this, where now both A and m a should be regarded as functions of time. Substituting the solution in Eq. 
( |2.3| ), we now obtain 

^-(m a A 2 ) = -3Hm a A 2 , (2.5) 

at 

or 

m a A 2 R 3 = constant, (2.6) 

where R is the scale factor of the universe, i.e., H = R/R. The energy density in these oscillations, at the 
present time, would be given by 

Po = \ (rn 2 a A 2 ) . (2.7) 

The axion mass in the present universe is determined by /pq. It depends mildly on the axion models, and 
is given by 

("ia)o/pQ - rn^U , (2.8) 

where the quantities on the right sides relate to properties of the pions, and are very well-known. However, 



we don't know the magnitude of the oscillations in the present universe. Therefore, we use Eq. (2.6) to write 

Po = l(m a ) (m a A 2 ) T x(R 3 /R 3 ), (2.9) 
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Figure 2: Summary of axion search results. The shaded regions are excluded. Taken from Ref. |2j. Some of 
the legends have been slightly edited to conform with the notations used in this paper. 



where R is the value of the scale factor at any arbitrary temperature T. The oscillations started, as we said, 



around T ~ T< 



QCD 



100 MeV. At that time, 



iJ ~ T /Mp, and A ~ /pq. Using these estimates and 



utilizing the fact that R/Rq = Tq/T, we obtain 



po^-- 



lmj, 



/PQ 



M P 



To V 
TqcdJ 



which gives 



0„ 



10 



-12 



/PQ 

lGeV 



(2.10) 



(2.11) 



Thus, axions can be an important source of dark matter if /pq is close to 10 12 GeV. 

Let us now see what is experimentally known about these axions. In Fig. @, we show the present status 
[[2|. The horizontal axis is the axion mass. The vertical axis denotes its coupling with two photons, which 
allows the decay mode a — > 77. As I said earlier, there is slight model dependence in these quantities. The 
predictions of two compelling models are shown in the graph by the lines marked DFSZ and KSVZ. 

If the axion mass is in the eV range, the decay photons could be detected through telescopes, and the 
failure to do so has been indicated in the figure. However, if /pq indeed equals 10 12 GeV and we take 
Eq. (pTq) as a strict equality, the axion mass comes out to be 10 -5 eV. So, if the axions has to contribute 
significantly to the energy density of the universe, we should find them around that value of mass. In this 
case, the two photon decay mode is so slow that the axion can be considered stable on cosmological time 
scales. However, because of the coupling, the axions may be resonantly converted into a microwave signal 
in a cavity where a strong magnetic field is present. These experiments are very promising. As we see, they 
are already probing very close to the interesting region. 



3 WIMPs 



We now come to the WIMPs, which is an acronym for "Weakly Interaction Massive Particles." There are 
a variety of ways in which such particles may be motivated. The most recent and most compelling of them 
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Figure 3: Summary of WIMP search results and future prospects. From Ref. 



is supersymmetry. Supersymmctric models have new particles — bosonic partners of known fermions and 
fcrmionic partners of known bosons. There are phenomenolgically strong reasons to suspect that these new 
particles, or "superpartners" , can be produced or annihilated only in pairs. If that is true, the lightest of 
these superpartners will be a stable particle. It will not decay spontaneously. These can be examples of 
WIMPs. There are other models also which predict WIMPs. 

The number density of these WIMPs in the early universe is governed by the Boltzmann equation. In 
the homogeneous and isotropic universe, this takes the form 



— = -3Hn - <rv(n 2 - n 2 ) . 
dt q 



(3.1) 



where n eq is the equilibrium number density for the temperature at any given time. The decrease of number 
density due to the expansion of the universe is represented in the first term in this equation, and the same 
due to pair annihilations goes in the second term. 

It is convenient to switch to dimensionless variables defined by 



n 

y= -■. 



m/T , 



(3.2) 



where m is the mass of the WIMP and s is the entropy of the background photons at a temperature T. The 
evolution equation now becomes 



®V _ ycq 1 ann 

dx x H(x) 



y 

Vcq 



1 



(3.3) 



where r aTln = 



aim ^cq^^ 



-3/2 



e x crv, putting in the equilibrium number density for a non-relativistic 



particle. The important thing to notice here is that y hardly changes if r an n "C H(x) at any given era 
characterized by a temperature determined by x. Using H(x) ~ T 2 /Mp = m 2 /Mpx 2 , we find 



H{x) 



mM P x 1 / 2 e x uv . 



(3.4) 



As x — * oo, i.e., for large time, this ratio goes to zero. It means that after a while, the pair annihilation rate 
will be negligible, because the universe has meanwhile become large enough so that the particles do not find 



each other. The remnant density of the WIMPs can be obtained by numerical integration of Eq. (3.3), but 
a rough estimate can be obtained by the density at the era when r ann = H(x). The solution of x Eq. (3.4) 
is now determined once the mass and the annihilation cross section are known. After this era, the number 
density just scales inversely as the volume of the universe. 

Let me now turn to the experimental searches j|. These are summarized it in Fig. |3[ As I just said, the 
parameters which determine the density of WIMPs in the present universe are the mass and the annihilation 
cross section. The two axes in this plot are precisely these two parameters. For the cross section, only 
the cross section with the nucleons is probed in the experiments. The solid lines are exclusions plots from 
already published data, and the dotted lines are projections for future experiments. The scatter plots are 
expected values of the parameters for various choices of the basic parameters in the supersymmetric model. 
It does seem that the experiments in the near future will be able to decide whether WIMPs really dominate 
the energy density of the universe. 
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